Optical memory based on ultrafast wavelength switching in a bistable microlaser 
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We propose an optical memory cell based on ultrafast wavelength switching in 
coupled-cavity microlasers, featuring bistability between modes separated by sev- 
eral nanometers. A numerical implementation is demonstrated by simulating a two- 
dimensional photonic crystal microlaser. Switching times of less than 10 ps, switch- 
ing energy around 15-30 fJ and on-off contrast of more than 40 dB are achieved. 
Theoretical guidelines for optimizing the performance of the memory cell in terms of 
switching time and energy are drawn. 
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In recent years, several designs of an optical memory (flip- 
flop) cell were proposed. Elements based on bistability in ac- 
tive laser-based systems iiL|2|,[3|,|4|] were shown to be promis- 
ing in terms of switching performance. All-optical switch- 
ing between the two memory states at the sub-nanosecond 
time- scale was reported for directionally bistable injection- 
locked microring lasers [2, 5] and for bistable gain-quenching 
semiconductor optical amplifiers [3]. Microlasers offer an ad- 
vantage in the switching energy due to a small cavity size 
and high Q-factors |6]. Recently, it was reported [H 0, [II] 
that a coupled-cavity microlaser can lase into either of its su- 
permodes in a bistable manner with a possibility of ultrafast 
picosecond-scale switching between these modes (80. 

In this Letter, we build up on our recent findings on bista- 
bility in microlasers |@, 0, 0] to propose a design of an ultra- 
fast optical memory cell. We consider two identical coupled 
defects in a two-dimensional (2D) triangular protonic crystal 
(PhC) lattice (see Fig. []} with active centers (e.g., quantum 
dots or quantum wires) embedded into the defects |H, [Toll . 
Such structures are within the state of the art of current fabri- 
cation possibilities ['11]. Direct SNOM observations confirm 
the formation of coupled-cavity supermodes in such a system 
fl2h . Unlike in polarization bistability, directional bistabil- 
ity in bidirectional ring lasers, or bistability among longitu- 
dinal modes in unidirectional ring lasers 11311 . bistable modes 
in coupled-cavity microlasers can be up to tens of nanometers 
apart from each other (8D, which makes the output signal sub- 
ject to an easy spectral separation compatible with a WDM 
arrangement. Our numerical simulations show that for a laser 
operating 100-300 times above threshold (which is still prac- 
tically realizable), switching between modes 3 nm apart can 
occur at less than 10 ps (see Fig. Q]), with switching energy 
estimated to be around 15-30 fJ. 

To describe the laser dynamics, we use the coupled mode 
theory [ 6], which is a generalization of existing models (see, 
e.g., [13]) for the case of arbitrary cavity shape. A res- 



onator comprising two identical coupled cavities supports two 
modes (the bonding and the antibonding supermodes) char- 
acterized by symmetric (in-phase) vs. antisymmetric (out- 
of-phase) cavity field distribution, respectively (see inset 
in Fig. [T]) We assume the supermodes to have frequen- 
cies C0\ 2 = <9o =F A© and decay rates K*i ; 2- The total cavity 
field E(r,t) can be written in terms of spatial mode distri- 
butions uj(r) and time-dependent mode amplitudes Ej(t) as 
E(r,t) = ui(v)Ei(t)er ia * + u 2 (r)E 2 (t)Q~ i0>2t . In the limit of 
class-B laser, the mode amplitudes obey the equations ((J 



-(Kj/2) (Ej+E] 



k 



(1) 



where /3i,2=A a ± iA w and £ « ji 2 cOo/2£ £h. The terms E? xt (t) 
account for the 7-th mode projection of the externally injected 
signal [14]. Njk(t) = £ j d 3 ru*j(r)N(r,t)uk(r) are projections 
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Fig. 1. (Color online) Analytical and numerical demonstra- 
tion of a 1 — » 2 —> 1 mode switching in a 2D PhC coupled- 
defect microlaser. E^ l (t) are switching signals in the defects; 
^1,2 (0 are FDTD output for the two modes; Ej(t) are ob- 
tained from Eqs. CQ)-0. 
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Table 1. The parameter values used in the calculations. 



Symbol Value Comment (see also [6] for details) 





~ l(r 28 A-s 


dipole moment of atomic laser transition 


£ 


12.3763 


dielectric constant of the host (GaAs) 


(Do 


2.13 xlO 15 s- 1 


operating frequency (around 885 nm) 


Aa 


1.88 x 10 13 s _1 


polarization decay rate (laser linewidth) 




0.1A a 


intermode frequency separation 


K j 


0.1 A a 


cavity mode decay rate 


Y 


0.001A a 


population inversion decay rate 


a jk 


3.51 x 10 11 


non-zero mode overlap integrals 



of population inversion Af (r, t) governed by the equations 

Njk = r(Rjk-Nj k )-^ ZM^N^Ej^-^'. (2) 

mnpq 

Here £ = ji 2 /4h 2 and Rjk are projections of the external 
pumping rate R(r,t) = R. Other parameters are listed in Ta- 
ble [T] The elements of the matrix M are composed of f3j and 
a mn _ £ j d 3 r^(r)^(r)w^(r)^(r). For symmetric pumping 
and near-orthogonal modes, the non-zero elements are 

M T PP = a T ( 1 /A» + 1 IK) , M 7' pq = 

M%>" m = (a]'f + a»])/f} m , M™' pp = a]F/p n . (3) 

It was shown that Eqs. ([T])-0 support bistable solutions pro- 
vided a J k l ~ oc^ and the pumping rate is high enough [6fl. A 
corresponding flow diagram is shown in Fig. [2^, clearly indi- 
cating two stable states (Ml and M2) with two equally sized 
domains of attraction divided by a separatrix. 

The principle of optical memory operation proposed here 
is based on the switching of the bistable microlaser between 
its two stable states by injection locking. Either mode can be 
locked into by choosing the right spatial profile of the injec- 
tion signal Ej xt while keeping the same spectral profile 
Suppose the laser is at M2, i.e., generating into mode 2 so that 
£2(0 = ^mode- An injected pulse shifts the mode balance in 
favor of mode 1, i.e., moves the resonator away from M2 and 
over the separatrix. Due to bistability, the system will relax 
to the other stable point Ml, i.e., ends up lasing into mode 1 
(see Fig. [2^). In order to reach the separatrix the minimal en- 
ergy delivered by the switching pulse W sw should be at least 
W c / \fl where W c is the energy stored in the cavity. 

In response to the switching pulse, assumed to be a step 
function with height E sw and duration cr, two processes can 
be distinguished (Fig.[2k): evolution of the system away from 
M2 under the action of the pulse and relaxation to Ml af- 
ter the pulse has subsided. The former is time-limited by the 
switching pulse duration a. The duration of the latter x r de- 
pends only on the initial point (h,h) on the flow diagram 
that the switching pulse brings the system to. This duration 
is large near the separatrix and rapidly decreases away from it 
(Fig.[2j)). It is thus desired that the switching pulse brings the 
system close to mode-1 axis or Ml. The desired pulse energy 
is then W sv/ ~ W C V2. When T¥ sw varies between W c /V2 and 
W c Vl, the switching time is expected to drop significantly. 



This can be confirmed by plotting the switching diagrams 
in (7i,/2)-space (Fig. [3]). One can see that while the pulse re- 
mains in effect, the system evolves along a straight line, build- 
up path (BP), accompanied by relaxation oscillations. The 
pulse duration a determines how far the system progresses 
along BP (I in Fig. [3^). As soon as the pulse ends, the system 
leaves BP at a drop-off point (DOP), quickly relaxing onto the 
M2-M1 line (II) and then slowly continuing along this line 
towards the other stable point (III). 

The slope angle (p of the BP depends on the balance be- 
tween £ sw and K. Higher £ sw increases the injection rate for 
the mode 1 while higher g-factors (smaller k*) make the de- 
cay of the mode 2 slower. Both effects lead to a decrease 
of (p (Fig. [3j)). Hence, setting the finesse of the resonator 
too high makes switching excessively long while lowering the 
g-factors increases R^ and, in turn, the operating pumping 
rate R. E sw and a can be selected so as to provide the fastest 
switching. Most of the BP should be covered and the DOP 
should be as close to Ml as possible (Fig. [3]). There is a clear 
tradeoff between switching time and energy, as seen in the 
plots T sw (£ sw , cr) in Fig. [3] A similar tradeoff was recently ob- 
served in low-energy all-optical switching between transverse 
patterns in atomic vapor and semiconductor microresonators 
El 0]. Further increase of E sw or o moves the DOP away 
from Ml, causing a slight growth of T sw . Operating the laser 
beyond the minimum in T sw may be useful to improve the re- 
liability if system parameters vary due to fabrication disorder 
or environmental factors. 

These optimization guidelines were used in the design 
of a proof-of-the principle memory cell for direct numer- 
ical simulations (Fig. Q]). The choice of geometry (cou- 
pled defects in a PhC lattice) was based on state-of-the- 
art fabrication capabilities full . The active medium was 
modeled using 4-level laser rate equations coupled to a 2D 
finite-difference time-domain (FDTD) solver [14]. Switch- 
ing is effectuated by in-phase (Ef xt = E^ xt ) and out-of-phase 
(Ef xt = —Ey Xt ) quasi-monochromatic rectangular pulses with 
frequency C0q delivered into the defects through PhC waveg- 
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Fig. 2. (Color online) (a) An above-threshold flow diagram for 
class-B two-mode microlaser described by Eqs. ([T])-© show- 
ing bistable lasing regime and the schematics of the mode 
switching; (b) density plot of time T r {hJi) it takes the sys- 
tem to relax from the point of origin {h,h) to the 20 dB- 
neighborhood of the corresponding stable point. 
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Fig. 3. (Color online) (7i,/2)-space diagrams of the 2 — > 1 
mode switching by a stepwise pulse: (a) with fixed E sw = 
0.5£ mo de and varying cr; (b) with fixed a = 3.5 ps and varying 
E sw . The insets show the dependencies t sw (£ sw ) and T sw (cr). 



uides. Due to symmetry matching, these patterns couple to 
the bonding (Ml) and the antibonding (M2) mode, respec- 
tively. The output radiation is extracted from another PhC 
waveguide (the "bus"). The signal for the two frequencies 
is determined as the across-the-bus flux of the Poynting vec- 
tor 3>i,2(0 ~ Jbus E ( r ' £ °i 5 2,0 x B*(r,fi)i j2 ,0 - nd2r obtained 
from coarse-graining Fourier transform of the fields at (Q\ 2 as 

E(r, ©1,2,0 = J^Eirj'y^t'dt' where T = 0.5 ps. This 
mimics the response of a spectrally selective time-resolving 
detector. The laser was pumped at 300 times threshold in or- 
der to assure bistable regime despite a considerable mismatch 
of the mode decay rates (K\ ^ K2), which is also the reason 
why #1 and <£ 2 show different on-off contrast 

Numerical simulations demonstrate switching between the 
two modes on the scale of several (less than 10) picoseconds 
(Fig. [U, with a decent qualitative agreement between FDTD 
and coupled-mode results. Analytical calculations tend to un- 
derestimate the switching times compared to numerics, both 
in this work and in earlier ones EH]- One possible reason 
may be due to the difference in laser system (4-level in FDTD 
vs. 2-level in analytics) but this remains a topic of further 
studies. The on-off contrast is at least 40 dB. To estimate the 
switching energy one should take into account that the energy 
transferred to a micro structure is a fraction p of the externally 
delivered power. In such lasers, the operating absorbed power 
P = 60 /iW and p is around 0.1-0.2. From the simulations, 
the cavity field enhancement factor was / = 10. Thus, for 
£" sw = £mode the estimation for the switching pulse energy is 



pfPo ~ 15-30 fJ, comparable to the earlier accounts fil. I t 
can be further reduced, e.g, using coupling optimization flivh . 
The storage time of the proposed device is theoretically infi- 
nite as long as the laser stays pumped and the in-cavity energy 
fluctuations remain smaller than W c / V2. 

The authors acknowledge helpful assistance of C. Kremers 
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the Deutsche Fouschungsgemeinschaft (FOR 557). 

References 

1. H. Kawaguchi, IEEE J. Sel. Top. Quant. Electron. 3, 
1254-1270(1997). 

2. Y. Liu, M. T. Hill, N. Calabretta, H. de Waardt, 

G. D. Khoe, and H. J. S. Dorren, IEEE Photon. 
Tech. Lett. 15, 1461-1463 (2003). 

3. J. Oksanen and J. Tulkki, Appl. Phys. Lett. 88, 181118 

(2006) . 

4. S. Osborne, K. Buckley, A. Amann, and S. O'Brien, 
Opt. Express 17, 6293-6300 (2009). 

5. M. T. Hill, H. J. S. Dorren, T. de Vries, X. J. M. Lei- 
jtens, J. Hendrik den Besten, B. Smalbrugge, Y.-S. Oei, 

H. Binsma, G.-D. Khoe, and M. K. Smit, Nature 432, 
206-209 (2004). 

6. S. V. Zhukovsky, D. N. Chigrin, and J. Kroha, 
Phys. Rev. A 79, 033803 (2009). 

7. S. V. Zhukovsky, D. N. Chigrin, A. V. Lavrinenko, and 
J. Kroha, Phys. Stat. Sol. B 244, 1211-1218 (2007). 

8. S. V. Zhukovsky, D. N. Chigrin, A. V. Lavrinenko, and 
J. Kroha, Phys. Rev. Lett. 99, 073902 (2007). 

9. K. A. Atlasov, K. F. Karlsson, E. Deichsel, A. Rudra, 
B. Dwir, and E. Kapon, Appl. Phys. Lett. 90, 153107 

(2007) . 

10. P. Gallo, M. Felici, B. Dwir, K. A. Atlasov, K. F. Karls- 
son, A. Rudra, A. Mohan, G. Biasio, L. Sorba, and 
E. Kapon, Appl. Phys. Lett. 92, 263101 (2008) 

11. K. A. Atlasov, K. F. Karlsson, A. Rudra, B. Dwir, and 
E. Kapon, Opt. Express 16, 16255-16264 (2008). 

12. S. Vignolini, F. Intonti, M. Zani, F. Riboli, 
D. S. Wiersma, L. H. Li, L. Balet, M. Francardi, A. Ger- 
ardino, A. Fiore, and M. Gurioli, Appl. Phys. Lett. 94, 
151103(2009). 

13. L. M. Narducci, J. R. Tredicce, L. A. Lugiato, 
N. B. Abraham, and D. K. Bandy, Phys. Rev. A 33, 
1842-1854(1986). 

14. S. V. Zhukovsky and D. N. Chigrin, Phys. Stat. Sol. B 
244, 3515-3527 (2007). 

15. R. Kheradmand, M. Sahrai, H. Tajalli, G. Tissoni, and 
L. A. Lugiato, Eur. Phys. J. D 47, 107-112 (2008). 

16. A. M. C. Dawes, D. J. Gauthier, S. Schumacher, 
N. H. Kwong, R. Binder, and A L. Smirl, Laser and Pho- 
ton. Rev., DOI 10. 1002/lpor.2008 10067 (2009). 

17. C. R. Otey, M. L. Povinelli, and S. Fan, Appl. Phys. Lett. 
94, 231109 (2009) 



